Search results for "Jacobian matrix"
showing 10 items of 47 documents
A new method for creating sparse design velocity fields
2006
We present a novel method for the computation of mesh node sensitivities with respect to the boundary node movement. The sensitivity field is sparse in a sense that movement of each boundary node affects only given amount of inner mesh nodes, which can result in considerable savings in the storage space. The method needs minimal control from the user, and it does not place any restrictions (such as block structure) on the mesh. Use of the method is demonstrated with a shape optimization problem using CAD-free parametrization. A solution to the classical die-swell free boundary problem by coupling the boundary node locations with the state variables is also presented. In that case, sparsity …
New Types of Jacobian-Free Approximate Riemann Solvers for Hyperbolic Systems
2017
We present recent advances in PVM (Polynomial Viscosity Matrix) methods based on internal approximations to the absolute value function. These solvers only require a bound on the maximum wave speed, so no spectral decomposition is needed. Moreover, they can be written in Jacobian-free form, in which only evaluations of the physical flux are used. This is particularly interesting when considering systems with complex Jacobians, as the relativistic magnetohydrodynamics (RMHD) equations. The proposed solvers have also been extended to the case of approximate DOT (Dumbser-Osher-Toro) methods, which can be regarded as simple and efficient approximations to the classical Osher-Solomon method. Som…
Locating Objects Away from Earth Surface: Positioning Accuracy
2013
The motion of the Galileo and GPS satellite constellations is simulated in Schwarzschild space-time, whereas photons travel in Minkowski space-time. This is a good enough approach to deal with the main goal of this paper: the study of positioning accuracy in the framework of the so-called relativistic positioning. Our study is based on numerical 4D simulations. In this meeting, the contribution of J.A. Morales-Lladosa contains some basic ideas which have been important to perform our numerical calculations. For four chosen emitters (satellites) of a certain constellation, many receivers located at different distances from Earth surface and in distinct directions are considered. Thus, we ver…
A Flux-Split Algorithm Applied to Relativistic Flows
1998
The equations of RFD can be written as a hyperbolic system of conservation laws by choosing an appropriate vector of unknowns. We give an explicit formulation of the full spectral decomposition of the Jacobian matrices associated with the fluxes in each spatial direction, which is the essential ingredient of the techniques we propose in this paper. These techniques are based on the recently derived flux formula of Marquina, a new way to compute the numerical flux at a cell interface which leads to a conservative, upwind numerical scheme. Using the spectral decompositions in a fundamental way, we construct high order versions of the basic first-order scheme described by R. Donat and A. Marqu…
Isotropic p-harmonic systems in 2D Jacobian estimates and univalent solutions
2016
The core result of this paper is an inequality (rather tricky) for the Jacobian determinant of solutions of nonlinear elliptic systems in the plane. The model case is the isotropic (rotationally invariant) p-harmonic system ...
Jacobian-Free Incomplete Riemann Solvers
2018
The purpose of this work is to present some recent developments about incomplete Riemann solvers for general hyperbolic systems. Polynomial Viscosity Matrix (PVM) methods based on internal approximations to the absolute value function are introduced, and they are compared with Chebyshev-based PVM solvers. These solvers only require a bound on the maximum wave speed, so no spectral decomposition is needed. Moreover, they can be written in Jacobian-free form, in which only evaluations of the physical flux are used. This is particularly interesting when considering systems for which the Jacobians involve complex expressions. Some numerical experiments involving the relativistic magnetohydrodyn…
Point counting on Picard curves in large characteristic
2005
We present an algorithm for computing the cardinality of the Jacobian of a random Picard curve over a finite field. If the underlying field is a prime field Fp, the algorithm has complexity O(p).
Spectral WENO schemes with Adaptive Mesh Refinement for models of polydisperse sedimentation
2012
The sedimentation of a polydisperse suspension with particles belonging to N size classes (species) can be described by a system of N nonlinear, strongly coupled scalar first-order conservation laws. Its solutions usually exhibit kinematic shocks separating areas of different composition. Based on the so-called secular equation [J. Anderson, Lin. Alg. Appl. 246, 49–70 (1996)], which provides access to the spectral decomposition of the Jacobian of the flux vector for this class of models, Burger et al. [J. Comput. Phys. 230, 2322–2344 (2011)] proposed a spectral weighted essentially non-oscillatory (WENO) scheme for the numerical solution of the model. It is demonstrated that the efficiency …
Mappings of finite distortion: Monotonicity and continuity
2001
We study mappings f = ( f1, ..., fn) : Ω → Rn in the Sobolev space W loc (Ω,R n), where Ω is a connected, open subset of Rn with n ≥ 2. Thus, for almost every x ∈ Ω, we can speak of the linear transformation D f(x) : Rn → Rn, called differential of f at x. Its norm is defined by |D f(x)| = sup{|D f(x)h| : h ∈ Sn−1}. We shall often identify D f(x) with its matrix, and denote by J(x, f ) = det D f(x) the Jacobian determinant. Thus, using the language of differential forms, we can write
A Fast Imaging Technique Applied to 2D Electrical Resistivity Data
2014
A new technique is proposed to process 2D apparent resistivity datasets, in order to obtain a fast and contrasted resistivity image, useful for a rapid data check in field or as a starting model to constrain the inversion procedure. In the past some modifications to the back-projection algorithm, as well as the use of filtering techniques for the sensitivity matrix were proposed. An implementation of this technique is proposed here, considering a two-step approach. Initially a damped least squares solution is obtained after a full matrix inversion of the linearized geoelectrical problem. Furthermore, on the basis of the results, a subsequent filtering algorithm is applied to the Jacobian ma…